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Abstract 



The two-point functions of the energy-momentum tensor and the Noether 
current are used to probe the 0(3) nonlinear sigma model in an energy range 
below 10 4 in units of the mass gap m. We argue that the form factor approach, 
with the form factor series trunctated at the 6-particle level, provides an almost 
exact solution of the model in this energy range. The onset of the (2-loop) 
perturbative regime is found to occur only at energies around 100m. 
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1. Introduction 



The 0(3) Nonlinear Sigma (NLS) model has long been appreciated as a 2-dimensional 
testing ground for nonabelian gauge theories. In particular it shares features like dynam- 
ical mass generation, running coupling constant and (conjectured) asymptotic freedom 
with its 4-dimensional cousins. A central issue is therefore to determine the domain of va- 
lidity of perturbation theory in the model. This requires an approach which is capable of 
describing the system both in the non-perturbative low energy and the perturbative high 
energy region. Although it could not have been expected from the outset, the form factor 
approach |], |2j (which essentially provides a low energy expansion of n-point functions 
based on the exact S"-matrix ||) turns out to be well suited for that purpose. Using the 
two-point functions of the energy-momentum tensor and the Noether current to probe the 
system, we find that the form factor series truncated at the 6-particle level provides an 
almost exact solution of the model in the energy range below 10 4 in units of the mass gap 
m. A similar fast convergence of the low energy expansion has previously been observed 
for the class of models with diagonal S'-matrices ||, §, |(| . 

We compare our results with Monte Carlo data [0 and 2-loop perturbation theory (PT). 
Within the range considered 2-loop PT appears to yield an accurate (within one percent) 
description of the system only for energies above 100m, provided one uses the known exact 
value of the Lambda parameter || to fix its absolute normalization. In 4-dimensional 
gauge theories of course neither the exact Lambda parameter nor the (almost) exact form 
factor curve is available. Based on the low energy Monte Carlo data alone one might 
therefore be tempted to maximize the apparent domain of validity of PT by tuning the 
Lambda parameter such as to match the relevant part of the Monte Carlo data. Doing 
this in the NLS model would give a value for the Lambda parameter that is smaller than 
the exact value by about 10%. This emphasizes the importance to have an independent 
estimate for the onset of the (2-loop) perturbative regime. 

2. Spectral representation of two point functions 

The form factors characterize an (integrable as well as non-integrable) QFT in a similar 
way as the n-point functions do. Assuming the existence of a resolution of the identity 
in terms of asymptotic multi-particle states, the n-point functions can in principle be 
recovered from the form factors. In the important case of the two-point functions this 
amounts to the well-known spectral representation. Explicitly for the Minkowski two- 



1 



point function (Wightman function) of some local operator O one obtains in a first step 

W°( x -y) = J2 — I IT — e~ i( - x °~ y ° )p o n)( - e) ~ i( - xl ' yl)p i n)( - e) \f ( °' n) (9)\ 2 , (1) 

>i iris ! J -i 4ti~ 

where /(°> n )(0) = (vac|0(O)|0 n , . . . , 6 X ) are the form factors of O and Pfr\6) = EiPM, 
Po(6) = mch.9, pi{9) = msh.9 are the eigenvalues of energy and momentum on an n-particle 
state.Q The local operators are classified by various quantum numbers, in particular by 
their Lorentz spin s and their mass dimension A. The form factors of a local operator 
with quantum numbers (A, s), \s\ < A can be assumed to be of the form 

A+s A — s 

f M W=i'(^T,e $ <y (^E^) 2 f in \0), (2) 

where f {n) {9) carries quantum numbers (A, s) = (0, 0) and is a function of the rapidity 
differences only. We shall henceforth always work with the 'scalarized' form factors f^ n \0) 
and drop the superscript l O\ For many purposes it is useful to rewrite (1) in the form of 
a Kallen-Lehmann spectral representation. Changing integration variables according to 

Ui = 9i — 9 i+ i , 1 < i < n — 1 , a = In 
M (n) (u) 2 = m 2 [n + 2j2 ch(u; + . . . + 

i<j 

and considering the case of an operator O with quantum numbers (A, s) one obtains 

(■OO 

W°(x-y) = -t dpp{p)D^ s \x-y-p), 
Jo 

\ji/m] 

M= E P W (/^), pV(p) = 6(p-m)\f^\ 2 , 

n=l 

p^\p) = jf ^ w f!r X \f (n) (u)\ 2 5(p - MW(ti)) , n>2. (3) 

Notice that no problem of convergence arises for the spectral density. First, each n- 
particle contribution exists because, for fixed p, the integrand has support only in a 

^^Our kinematical conventions are: (x ^ 1 ) are coordinates on 2-dimensional Minkowski space M 1 ' 1 
with bilinear form x ■ y — x p "r\ lLV y v , r) = diag(l, — 1). Lightcone coordinates are — (x° ix 1 )/^/^ = x^. 
The normalization of the 1-particle states is (0i\02) = 4tt6(0i — 02 ), which corresponds to the standard 
normalization in d + 1 dimensions, specialized to d = 1. For simplicity we suppress internal indices for 
most of this section. 
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compact domain V((J,) C IR™ -1 in which the form factors are bounded functions, so that 
the integration is well-defined. Viewed as a function of // one observes that p^{jj) has 
support only for fi > mn. Therefore only a finite number of terms (those with n < [fj,/m], 
[x] being the integer part of x) contribute to p(/i). Under some mild assumptions on the 
growth of the existence of the spectral density then guarantees that of the 2-point 
function, as defined through (3). The integration kernel D^ A ' s \x;m) in (3) is given by 
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J 2,7V 

= (td + ) A+s (td^) A - s D(x; m) , (4) 

where D(x;m) := D^ '°\x; m) is the two-point function of a free scalar field of mass m. 
(In general however D^ A,s \x;m) does not coincide with the two-point function of some 
composite operator O built from the scalar field). In the spinless case one obtains from 
the wave equation (x; m) = (m 2 /2) A D(x; m). The support properties of the kernel 

functions can be extracted by evaluating them in terms of Bessel functions; in particular 
for strictly spacelike distances one finds 

(a; m) = i (^pj -L (^pj K 2s (mV^) , x 2 < . (5) 

Essentially this also yields the kernel function entering the spectral representation of 
the Euclidean two-point function (Schwinger function). The latter can be defined by 
S°(xi,x 2 ) := W°(—ix2,Xi) for x 2 > and then by analytic continuation to x 2 < 0. In 
the spectral representation (3) this yields —iD^ A ' s \—ix 2 ,xi,m) for the Euclidean kernel 
function with the right hand of (5) taken for D^ A ' s \x; to). Notice that the spectral density 
is the same in the Minkowski and in the Euclidean case. Return now to the Minkowski 
space situation. For comparison with perturbation theory one needs the time-ordered 
two-point function and its Fourier transform. Its spectral representation is easily read off 
from (3), (4) 

G°{x -y) = -i dp p(/i) G (A > s \x -y-ti), 
Jo 

G {A ' s \x; to) = 9(x ) D^ s \x; to) - 9(-x ) [D^ s \x; to)]* 

f _fP_ p -i P - x (P + ) A+S (P-) A ~ S (6) 
J (2vr) 2 p 2 -m 2 + ie 1 ' 

For the Fourier transform this means 
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The definition of I(z) was chosen such that it has a cut along the negative real axis and 
one can recover the spectral density from the discontinuity along this cut. 

As an example consider the case of the energy-momentum (EM) tensor. It is of additional 
interest, because its spectral density is closely related to the Zamolodchikov C-function 
0. The spectral representation of the Minkowski 2-point function is 



poo 

w ™a( x -y) = -i J o dp pO) D mpa {x -y;ti), 

— 9(p )5(p 2 - m 2 ) e~ ip ' x (ri^p 2 - PpPv)(r] pcT p 2 - P P Po 



The lightcone components correspond to the SO (I, 1) -irreducible pieces; in particular 
D ++!++ (x; m) and _D + _ i+ _(:e; m) coincide with D^ 2,2 \x; m) and D^ 2 ^{x; m), respectively. 
Combining (3) and (5) one finds for the behaviour at small spacelike distances 

2 



X 



X 



+ 



+ 0' 1 



27T 2 (x 2 ) 2 \X 2 



x 2 < 



roc 

YItx \ dp pip) . 
Jo 



The number c as defined here coincides with the central charge of the Virasoro algebra 
in the conformal field theory describing the UV fixed point of the renormalization group. 
This latter fact is part of the statement of Zamolodchikov's C-theorem.Q The proof of the 
C-theorem is particularly transparent from the viewpoint of the spectral representation 
10 1 . The normalization is such that H = Jf^ dx l Tqo(x) has eigenvalue \/p 2 + m 2 on an 



asymptotic single particle state, where m is the mass gap. For the 2-particle form factor 
of the T + _ component this implies 

f( T +-> 2 \9 + in,9) =m 2 , (9) 

which fixes the normalization of all higher particle EM form factors. The scalarized form 
factors of all components can be interpreted as the form factors of a single field r defined 
by Tp U = E Pj pt va d p d a T. The conservation equation is then built into the parametrization 
(2). 

Also a Noether current comes with an intrinsic normalization, which arises from the Lie 
algebra of its conserved charges. In the case of the 0(3) NLS model one has [Q a , Qb] = 
itabc Qc and 

Qa \Q)b = itabc \0) c , 

2 In particular c is finite; for other operators the integral over the spectral density will in general not 
converge. 
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where Q a = J dx 1 Jo ia (x) is the conserved charge of the 0(3) Noether current. For the 
2-particle form factor of J$ A this implies 

f£- a (9 + in,9) = -2ime abc c\i6 , (10) 

which fixes the normalization of all higher current form factors. The scalarized form 
factors of both (spacetime) components can be interpreted as the form factors of the field 
r a defined through J^ a = d v r a . The two point function of J^ a (x) and J u ,b(y) and 
hence the spectral density is of course proportional to 8 ab . 

3. Polynomial 0(3)-irreducible form factors 

Solutions of the form factor equations with the 0(3)-invariant S'-matrix || have been 
obtained by Kirillov and Smirnov JTTI] by a fusing procedure from the previously known 
form factors of the SU(2) Thirring model ||. Unfortunately these solutions refer to a basis 
in the tensor product 3® n (21 + 1, 1 6 IN being the 0(3) irreps of spin I) that is related to 
the canonical basis through a complicated rapidity-dependent basis transformation. Since 
the most interesting local operators in the model all transform irreducibly under 0(3) one 
has to determine the intertwiners 

g(n) . g^n > 2 j + j 

to the canonical basis in (D 2 ' +1 = 21+1. If one starts from the Kirillov- Smirnov basis 
in 3® n also these intertwiners will be rapidity-dependent and their explicit computation 
is almost as difficult as the computation of the form factors themselves. We therefore 
found it easier to start afresh and to work with the canonical basis in 3® n and constant 
intertwiners. A detailed exposition of the technique and the results will be given elsewhere 

Ell- 
in survey we studied the following four local operators, which seem to be the most in- 
teresting ones in the 0(3) model: The 'fundamental' spin field S, the Noether current, 
the EM tensor and the topological charge (TC) density. For orientation we tabulate the 
expressions in terms of the field S, their quantum numbers (/, A) and the particle numbers 
of the form factors computed. 
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0[S] 

L J 


(I, A) 


ff's computed 


spin 


S 


(1,0) 


1,3,(5) 


current 


j fM = S x d^S 


(1,1) 


2,4,(6) 


EM tensor 


T liV = d li S-d u S-\<Or,di>S-d p S 


(0,2) 


2,4,(6) 


TC density 


q = S ■ (d^S x d v S) 


(0,2) 


3,(5) 



The explicit formulae can be found in [|l!J. The brackets around the last particle number 
indicate that the resulting expressions for the form factors are too long to be communi- 
cated in print; instead we shall give the expressions for their modulus squares, which is 
also the quantitiy entering the spectral densities.^ Higher particle contributions could be 
computed in principle, although more refined techniques would be necessarry to circum- 
vent limitations in computer power. It seems however that higher particle contributions 
to the spectral density are strongly suppressed in the 0(3)-model. Indeed, in section 4 
we shall argue that truncating the form factor series at the 5 or 6 particle level provides 
an 'almost exact' solution of the model up to energies of about 10 4 times the mass gap. 
'Almost exact' means that within this energy range the deviation of the true spectral 
densities from the truncated ones is estimated to be less than a percent. If one accepts 
the conventional wisdom about the asymptotic freedom of the model at very high ener- 
gies, the combination of the non-perturbative form factor results (for energies up to 10 4 
in units of the mass gap) and the perturbative expansion at energies above 10 4 would 
provide a solution of the model sufficient for most practical purposes. For issues of prin- 
ciple it would nevertheless be highly desirable to get some control on the higher particle 
contributions, in particular because the 'conventional wisdom' has been challenged fl3| . 



Some first results in this direction will also be presented in [12 . An important such issue 
is trying to actually prove (or disprove) asymptotic freedom. 

We shall not discuss the construction of the form factors here. Let us just mention 
the decomposition into 0(3) irreducible components. Let fi,A(@) denote the component 
of /a(0) in 21 + 1. Explicitly fi,A{&) is an n-index 0(3) tensor parametrized by mi(n) 
independent functions, where mi(n) is the multiplicity with which 21 + 1 occurs in 3® n . 
The choice of a parametrization amounts to a choice of an intertwiner Q\ between the 
0(3) representations 3® n and 21+1. Picking a basis in both spaces one can write 



fi,A(o) = Q? A g a (e) . 



3 The difference in size is enormous. The 6-particle form factor of the current has about 3.4 Mbytes, 
i.e. approx. 700 A4 pages; the most compact form of the square requires less than one page. 
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As remarked before, in contrast to iTTll, we work with the canonical basis in 3® n and 



constant (rapidity-independent) intertwiners Q^a- (This means that each component Q?a 
can explicitly be written as a combination of antisymmetric epsilon tensors and Kronecker 
deltas.) For the components g a {9) we used the following Ansatz 

9a(0) = g a (9) n 3 ^ 2 ]J ij(9 k - 9 r ) , 1 < a < m { (n) , 

fc>r 

The main structural result can then be summarized as follows [O, Il5l |12|: 



Fact: The form factor equations for the 0(3) NLS model decompose into decoupled recur- 
sive systems for the 0(3)-irreducible components g^ n \9) = (gi, . . . , g m ^)(9). There exist 
unique sequences of solutions (g^ n \0)) n >i of the form (12) s.t. g ( - n \9) = (g 1 , . . . ,g m ^ n )){9) 
are polynomials in the rapidities of total degree N = \{n 2 — 3n) + iVo and partial degree 
p = n - 2 + p . 

Here 'partial degree' means the degree in an individual variable. For the operators consid- 
ered before the initial values are given by (N ,p Q ) = (0, —1) for EM tensor & TC density 
and (iVo,po) = (1, 0) for Current & Spin. The fact that the form factors can be reduced 
to polynomial expressions is specific for the 0(3) model and will cease to hold for the 
O(N) NLS models with N > 3. The first part of the above statement in particular means 
that the 0(3) spin is a good quantum number for form factor sequences: If one member 
of the sequence is 0(3)-irreducible of spin / then all other members will be too. This is 
only seemingly trivial; it will not hold for any Lie group and any S'-matrix. 

For the computation of the spectral densities and the two point functions the modulus 
square of the form factors is needed. A version of \fi n \0)\ 2 that allows analytic continu- 
ation to complex rapidities is 

fl,a n ...ai(9n, ■ ■ ■ , @l) fl,ai...a n (#1, • • • , #n) = 9a{^ni ■ ■ ■ > #l) Ca,/9 9/3[dl, ■ ■ ■ , 9 n ) 

= a ivr 3 - 2 cf V) n m - Br) n m - ^ , (is) 

k>r i<j 

where C a ^ is a constant 'metric' and I = 0, 1 correspond to the EM tensor & TC density 
and Current & Spin series, respectively. The prefactor Ci is given by Co = 3 and C\ = 5 a b- 
The so defined square is a completely symmetric function of the rapidities; in particular 
its polynomial part G l (9) is a completely symmetric, boost invariant polynomial in the 
rapidities. (It depends only on the rapidity differences.) By choosing an appropriate basis 
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in the space of symmetric polynomials it can be brought into a very compact form, as 
compared with what one might guess from the size of the form factors (c.f. footnote 3). A 
convenient basis is obtained as follows: For fixed n let <j[ n \ . . . , denote the elementary 
symmetric polynomials in 9i, . . . , 9 n , i.e. 



a 



(n) _ 

k Z_^i u ii ■ ■ ■ u ik 

h<...<i k 



Oh ■ ■ ■ @i, 



Then define polynomials rj^ by 



rj: n \9 1 ,...,9 n ) = 4 n \9 1 ,...,9 n ) , 2 < k <n , 

r[ n) = -at ] , 9 3 = 9 j - -(9, + ... + n ) . (14) 
n n 

All monomials in the t^, k > 2 are manifestly boost invariant. (The price to pay is 
that the partial degree is no longer manifest.) For the sake of illustration we display the 
polynomial part of the 4-particle squares of the EM tensor and the Current in this basis 
(where the superscripts '(4)' are suppressed): 

EM tensor: 

G ( o\9) = 4(r 2 2 + 12r 4 ) - 7r 2 32r 2 + 28vr 4 , 
Current: 

gJ 4 V) = -4(6t 2 3 + 9r 2 + 40r 2 r 4 ) + 87r 2 (25r 2 2 + 44r 4 ) - 4487r 4 r 2 + 272tt 6 . (15) 



4. Results for spectral densities and two-point functions 

Knowing the form factor squares the evaluation of the spectral densities and the two-point 
functions is in principle straightforward. The integrations in (3) can be done numerically 
to good accuracy. Here we shall restrict attention to the EM tensor and the Noether 
current. For the spectral densities an accuracy of 10~ 3 was used and the results for the 
EM tensor and the Noether current are shown in figures 1 and 2, respectively. In the 
EM case one can use equation (8) to compute the n-particle contributions to the central 
charge. One finds 

POO 

c (2) = 1.603, c (4) = 0.194, c (6) = 0.072, c (n) = 12vr / dfxp {n \fi). (16) 

Jo 
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If one accepts asymptotic freedom as a working hypthesis, (tree level) PT predicts c = 2, 
corresponding to the two unconstrained bosonic degrees of freedom. The form factor 
computation shows that this is compatible with the non-perturbative low energy dynamics 
of the model and that even for an extreme UV quantity like the central charge the low 
particle contributions dominate. An alternative non-perturbative consistency check is 
provided by the thermodynamic Bethe Ansatz, which also yields c = 2 [16| . Let us point 
out that the value c = 2 does not carry much information about the nature of the UV 
limiting CFT. For a number of reasons it cannot be that of two decoupled free bosons. 



rho [x] 



0.012 



. 01 



. 008 



. 006 



. 004 



0.002 




Figure 1: 2-, 4-, and 6-particle contribution to the spectral density of the EM tensor as a 
function of x = log 2 (/z/m). 



Figures 1 and 2 illustrate some general features of the n-particle contributions to the 
spectral densities. First one observes a strong suppression of the higher particle con- 
tributions, as far as the value of the maximum and the enclosed circci Eire concerned. 
Generally speaking the maximum of p^ n \p) is smaller by 1.5 to 2.5 orders of magnitudes 
compared to the maximum of p("~ 2 \p), while the position of the maximum is shifted 
towards higher energies. Nevertheless at sufficiently high energies the n-particle contribu- 
tion overtakes the (n — 2)-particle contribution, i.e. p^ n \p) > p^ n ~ 2 \p), for p > p n , where 
P {Pn) = P^ n ~ 2 \pn)- (This feature is not directly visible in figures 1 and 2 because of 
the relative rescaling and the insufficient magnification of the intersection region.) The 
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point of intersection p n is of particular interest because it provides an intrinsic mea- 
sure for the quality of the approximation made by truncating the form factor series at 
the n-particle term: Since p n +2 the (n + 2)-particle contribution can safely be 

ignored up to energies p ~ p n in the sense that the correction to p^ (//) + ... + p^ (p) 
would be at most (of the order of) a percent. The results for these points of intersection 
(/i„,p (n_2) (/i n ) = p {n) (p n )) are as follows: 

EM: n = 4: (1.7 • 10 2 , 4.7 • 10" 6 ) , n = 6: (4.5 • 10 5 , 2.5 • 10" 10 ) , 
Current: n = 4: (1.6 • 10 2 , 3.8 • 10~ 4 ) , n = 6: (1.2 • 10 6 , 3.2 • 10~ 8 ) . (17) 

Thus, the form factor series truncated at 6 particles should provide accurate results for 
the spectral density up to energies O(10 5 ) - results, which can then be compared with 
those obtained by other techniques. 

In the high energy regime one expects the onset of asymptotic freedom and the pertur- 
bative predictions should coincide fairly well with the form factor curves. This is indeed 
the case as can be seen in figure 3. What is shown is the Fourier transform I(p 2 ) of the 
2-point function of the Noether current, computed once in 2-loop PT and once via (7) 
by truncation of the form factor series. In performing the integral transformation (7) one 



10 



can see that accurate results for p(/x) in the region 2 ~ [i ~ 10 5 produce accurate results 
for in the range z ~ 10 4 . In the energy range 100 ~ fx/m ~ 10 4 the coincidence of 
the 2-loop PT curve and the 2 + 4 + 6-particle form factor curve indicates that the system 
is desribed well by 2-loop PT in this regime. The result (17) for the intersection point hq 
suggests that the deviation at yet larger energies can entirely be attributed to the trunca- 
tion of the form factor series. Notice however that the onset of the (2-loop) perturbative 
regime occurs at much higher energies fi/m ~ 100 than is sometimes pretended in the 
4-dimensional counterpart of this situation. Figure 4 is a magnification of the low energy 
region of figure 3, where non-perturbative effects are expected to become important. 
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Figure 3: Comparison: Form factor approach versus 2-loop perturbation theory; logplot 
of I(p 2 ) against p/m. 



Here Monte Carlo simulations provide an alternative non-perturbative technique to probe 
the system [17 , |7|. The simulations were made using a Wolff-type cluster algorithm 
on a 460 square lattice at inverse coupling (3 = 1.80 (correlation length £ = 65.05). The 
agreement between the MC data and the form factor results is excellent. One also sees 
that for energies between 30 and 45 the PT curve runs almost parallel to the MC data. 
Without the guidance of the form factor result one would thus be tempted to match 
both curves by tuning the Lambda parameter approriately. Doing this however, the 
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Lambda parameter comes out wrong by about 10% (from below), as compared with the 
known exact result |§. Generally speaking one sees that a determination of the Lambda 
parameter from MC data and (2-loop) PT about the onset of the (2-loop) perturbative 
regime enters. The form factor results in the 0(3)-NLS model show that this regime sets 
in only at relatively high energies /i/m ~ 100, which should at least be taken as a warning 
in 4-dimensional gauge theories. 
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